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LONG-TERM  GOALS 

Develop  elastic  parabolic  equation  (PE)  method  capabilities  in  order  to  study  elastic  propagation 
mechanisms  and  their  effects  on  underwater  acoustic  environments  in  the  form  of  scattering  at  an  elastic 
interface,  oceanic  T -waves,  and  Scholte  waves. 

OBJECTIVES 

To  implement  explosive  and  earthquake  type  seismic  sources  in  current  elastic  parabolic  equation 
methods  for  underwater  acoustic  environments.  These  will  be  used  to  investigate  propagation 
mechanisms  that  depend  on  ocean  bottom  elastic  properties  and  result  in  acoustic  propagation  in  the 
ocean.  Two  examples  of  such  propagation  are  oceanic  T -waves,  which  are  acoustic  waves  that  result 
from  earthquake  or  buried  explosive  sources,  and  Rayleigh-type  waves  along  the  ocean  floor,  whose 
existence  requires  horizontal  and  vertical  displacements  present  in  elastic  sediments.  The  generation, 
propagation,  and  potential  contribution  of  these  waves  to  the  otherwise  quiet  acoustic  field  of  the  deep 
ocean  all  require  study,  in  particular  as  potential  sources  of  unexplained  deep  shadow  zone  arrivals  that 
have  been  experimentally  observed  below  the  ray-theoretic  turning  point.  [1,  2] 

APPROACH 

In  a  cylindrically  symmetric  environment,  where  r  is  the  distance  from  the  source  and  z  is  depth,  recent 
parabolic  equation  methods  for  acoustic  propagation  in  elastic  sediments  are  based  on  the  (ur,w) 
formulation  of  elasticity,  where  ur  is  the  horizontal  derivative  of  the  horizontal  displacement  and  w  is 
the  vertical  displacement.  [3]  The  outgoing  portion  of  the  separated  Helmholtz  operator  leads  to  the 
parabolic  equation  for  a  range-independent  environment, 


du 


(1) 


where  L  and  M  are  matrices  containing  depth-dependent  operators  that  incorporate  compressional  wave 
speed,  cp,  shear  wave  speed,  cs,  and  density  p  via  the  Lame  parameters  of  the  elastic  medium,  X  and  fi. 
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Range  dependence  is  included  by  modeling  sloping  boundaries  with  a  stair-step  approximation  and 
applying  appropriate  matching  conditions  at  each  vertical  interface.  This  approach  is  combined  with 
rotated  variable  methods[4]  and  an  elastic  single-scattering  approximation [5]  to  obtain  computationally 
stable  solutions  with  improved  accuracy  in  complex  multi-layered  range-dependent  underwater  acoustic 
environments. [6]  Parabolic  equation  self-starters  for  sources  in  a  fluid[7]  or  in  an  elastic  material[3] 
have  recently  been  tested  for  accuracy  in  underwater  acoustic  environments.  [8]  These  are  used  for 
analysis  of  the  effects  of  elastic  parameters  and  underwater  bathymetry  on  acoustic  and  seismoacoustic 
propagation. 

Wavenumber  content  of  range  independent  portions  of  an  elastic  PE  solution  f(r)  is  calculated  using 
the  Hankel  transform 

F(k,z)  =  f  f(r,z)Jo(kr)rdr  (2) 

Jr0 

where  ro  and  R  mark  the  beginning  and  end  of  a  range-independent  portion  of  the  environment,  and 
Jo(kr )  represents  the  Bessel  function  of  the  first  kind  of  order  0.  These  results  are  used  to  resolve  the 
nature  of  propagating  acoustic  modes  and  potential  interface  wave  modes. 

Downslope  conversion,  which  occurs  when  a  sloping  ocean  bottom  causes  the  resulting  transmission 
angle  of  incident  elastic  waves  to  focus  acoustic  waves  into  the  SOFAR  channel,  is  one  of  two  primary 
mechanisms  of  oceanic  T -wave  generation.  [9].  The  second  mechanism  generates  so-called  abyssal 
T -waves  from  areas  of  the  the  deep  ocean  where  there  is  no  significant  sloping  bottom.  It  is  believed 
that  ocean  bottom  roughness  scatters  the  elastic  waves  up  into  the  water  column  and  subsequently  into 
the  SOFAR  channel. [10]  Characteristics  of  underwater  acoustic  fields  in  the  presence  of  Normally 
distributed  ocean  bottom  roughness  are  determined  by  analysis  of  transmission  loss  results  and 
wavenumber  spectra  resulting  from  multiple  realizations.  [1 1] 

Scholte  interface  waves  are  excited  by  seismic  sources  and  have  been  observed  by  seismometers  at  the 
ocean  bottom.  [12,  13]  Energy  from  interface  waves  has  also  been  detected  by  hydrophones  near  the 
seafloor[13]  well  below  the  SOFAR  channel,  suggesting  these  waves  could  influence  deep-shadow  zone 
arrivals  observed  during  NPAL.[1,  2]  Elastic  wave  theory  predicts  the  Scholte  wave  speed  to  be 
approximately  0.8c5.  Thus,  a  peak  in  the  wavenumber  spectra  calculated  from  elastic  PE  solutions  near 
this  speed  indicates  interface  waves  have  been  excited. 

WORK  COMPLETED 

•  The  capability  of  PE  solutions  to  generate  abyssal  oceanic  T -waves  from  seismic  sources  has 
been  verified  by  comparing  transmission  loss  results  for  a  flat  seafloor  to  those  from  multiple 
Normally  distributed  rough  bottom  realizations. [14]  These  results  have  been  utilized  to  describe 
effects  of  rough  bottom  ocean  depth  on  T  wave  generation  and  propagation  by  examining 
depth-averaged  acoustic  transmission  results  from  elastic  PE  solutions.[ll] 

•  Generation  and  SOFAR  channel  propagation  of  oceanic  T  waves  from  seismic  sources  in  the 
presence  of  intervening  seamounts  or  coral  reef  barriers  is  established  using  elastic  PE 
solutions.[ll] 

•  Hankel  transforms  characterize  how  the  range  of  large  ocean  bottom  topography  affects  the 
presence  and  amplitude  of  Scholte  wave  arrivals  at  receiver  distances  of  100  km.  [14] 
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•  Wavenumber  spectra  from  a  one  layered  elastic  bottom  environment  with  an  intervening 
seamount  demonstrate  effects  of  elastic  parameters  on  relationship  between  water  column 
acoustic  energy  and  ocean  bottom  elastic  energy,  including  interface  wave  propagation. [15,  1 1] 

•  Elastic  PE  solutions  with  seismic  sources  [8]  determine  the  connection  between  slope  angle  and 
amplitude  of  propagating  oceanic  T -waves.  [1 1] 

•  Elastic  PE  solutions  for  long  range  SOFAR  channel  propagation  of  abyssal  T -waves  in  a  layered 
continental  shelf  environment  exhibit  Scholte  interface  waves  and  illustrate  that  these  solutions 
provide  on-shore  seismic  receptions.  [1 1] 

RESULTS 

The  generation  of  abyssal  oceanic  T -waves  by  small-scale  range  dependence  was  established  by 
demonstrating  acoustic  transmission  loss  averaged  over  several  realizations  of  Normally  distributed 
bottom  roughness  is  significantly  lower  than  acoustic  transmission  loss  for  perfectly  flat  ocean  bottoms. 
The  left  panel  of  Figure  1  shows  elastic  PE  results  for  a  4  km  deep  water  layer  with  a  Munk  deep  water 
sound  speed  profile  lies  on  top  of  an  elastic  halfspace  A  10  Hz  seismic  source  is  placed  2  km  below  the 
ocean  bottom  interface.  The  solid  red  curve  shows  transmission  loss  at  zr  =  1300  m,  which  is  the  axis 
of  the  SOFAR  channel.  The  solid  black  curve  is  the  average  transmission  loss  for  10  realizations  of 
Gaussian  ocean  bottom  roughness,  which  occurs  for  the  first  25  km  of  propagation.  The  shaded  region 
represents  one  standard  deviation  of  transmission  loss  results.  Until  about  40  km,  the  results  are  similar. 
However  as  range  increases,  additional  acoustic  energy  scatters  into  the  SOFAR  channel  and  begins  to 
dominate  the  propagation.  At  120  km  the  average  of  the  transmission  loss  from  the  rough  bottoms  is 
more  than  20  dB  for  both  compressional  and  shearing  seismic  source. 


Figure  1:  Elastic  parabolic  equation  results  for  4  km  deep  water  layer  over  an  elastic  half  space  with 
a  10  Hz  source  is  at  zs=6  km.  (Left)  Red  curves  show  transmission  loss  at  zr=  1300  m  when  the  ocean 
bottom  is  flat.  Black  solid  curves  shows  show  average  transmission  loss  at  zr=1300  m  for  10 
realizations  where  the  first  25  km  of  the  ocean  bottom  is  rough.  Shaded  region  shows  standard 
deviation  of  results  from  these  realizations  and  confirms  that  roughness  near  the  seismic  source 
leads  to  reduced  transmission  loss  associated  with  T -waves.  (Right)  Wavenumber  spectra  for  single 
rough  bottom  realization  in  the  same  environment  for  every  200  m  depth  show  evidence  of 
propagating  modes  in  the  water  column.  Depth  average  (red  curve)  and  standard  deviation  (shaded 
region)  of  wavenumber  spectrum  are  shown  at  bottom  of  panel. 
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The  right  panel  of  Figure  1  shows  wavenumber  spectra  at  120  km  range  and  several  depths  for  one 
realization  in  the  environment  described  above.  The  arrival  pattern,  specifically  the  peaks  near  the  blue 
dashed  line  (indicating  average  ocean  sound  speed)  to  the  right  of  the  figure  at  1250  m,  confirms  the 
presence  of  low-order  modes  propagating  in  the  SOFAR  channel.  Many  peaks  appear  at  lower 
wavenumber  values  (to  the  left  of  the  figure)  indicating  that  high  order  modes  remain  in  the  received 
signal.  These  peaks  are  particularly  strong  near  the  ocean  bottom  and  their  associated  signals  would 
arrive  much  later  than  those  observed  near  the  SOFAR  channel  axis. 


Figure  2:  (Left)  Elastic  parabolic  equation  solution  for  a  10  km  deep  10  Hz  shear  source  showing 
abyssal  T -wave  energy  traveling  in  the  SOFAR  channel  and  upslope  onto  land  via  a  beach.  The 
solution  includes  10  km  of  zero  water  depth  beginning  at  r=190  km.  Three-layer  ocean  bottom 
consists  of  elastic  sediment  layer,  a  transitional  elastic  layer,  and  a  hard  elastic  layer.  Normally 
distributed  roughness  occurs  along  the  first  25  km  of  the  ocean-bottom  interface.  There  is  an 
apparent  interface  wave  between  elastic  layers  (approx.  4250  m  depth )  and  evidence  of  mode  cutoff 
as  the  acoustic  energy  moves  upslope.  (Right)  A  5  Hz  shear  self-starter  is  20  km  below  the  ocean 
bottom  in  an  environment  featuring  a  seamount.  A  T  -wave  is  generated  by  the  downward  slope  of  the 
seamount  and  travels  to  a  range  of 400  km  in  the  SOFAR  channel  of  a  Munk  deep  ocean 

sound-speed  profile. 


The  left  panel  of  Fig.  2  shows  an  elastic  PE  solution  for  a  10  km  deep  10  Hz  seismic  source  in  a  4  km 
deep  ocean  environment  with  a  Munk  deep  water  sound  speed  profile.  Ocean  bottom  roughness  occurs 
on  a  layered  elastic  bottom  for  the  first  25  km  of  acoustic  propagation.  The  oceanic  T -wave  is  clearly 
propagating  and  has  multiple  turning  points  before  the  environment  begins  to  slope  upwards  as  for  a 
continental  shelf.  This  solution  illustrates  the  ability  of  the  elastic  PE  to  generate  abyssal  T -waves,  and 
apparent  interface  waves,  which  can  be  seen  along  the  interface  between  two  elastic  layers  at  4250  m 
depth.  Propagation  of  the  T -wave  as  purely  elastic  energy  can  be  seen  10  km  past  where  the  water  layer 
disappears. 

The  right  panel  of  Fig.  2  shows  a  solution  for  a  5  Hz  seismic  self-starter  20  km  below  the  ocean  bottom 
to  represent  a  very  deep  earthquake-type  source.  The  ocean  layer  is  2.5  km  deep  at  r  =  0  and  lies  over  a 
layered  elastic  ocean  bottom.  The  top  of  the  seamount  intersects  the  ocean  surface,  assuring  that 
acoustic  energy  to  the  right  of  the  seamount  has  been  converted  from  elastic  waves  in  the  seamount  via 
downslope  conversion.  Substantial  acoustic  energy  is  transferred  to  the  oceanic  T -wave,  which 
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propagates  near  the  SOFAR  channel  axis  (approx.  700  m  deep  in  this  example)  to  a  range  of  400  km. 
This  solution  could  be  extended  to  much  longer  ranges. 


Figure  3:  (Left)  Effects  of  changing  slope  angle  on  depth  averaged  transmission  loss  at  150  km 
range  for  a  5  Hz  compressional  (solid  line  with  circles)  and  shear  (dashed  line  with  squares)  seismic 
source  at  10  km  depth.  The  angle  of  0°  represents  a  4  km  deep  flat  environment  In  both  cases  the 
transmission  loss  values  show  a  trend  that  decreases  with  increasing  slope  angle.  (Right) 
Wavenumber  spectra  at  zr=2450  m  for  varying  cp  in  a  2500  m  deep  ocean  environment  with  a 
one-layered  elastic  bottom  cs  =  0.5 cp.  A  10  Hz  acoustic  source  is  located  at  zs=500  m.  A  2000  m  tall 
seamount  lies  between  the  source  and  the  receiver,  which  is  100  km  away.  Large  amplitude  of 
propagating  modes  (near  blue  dashes)  corresponds  to  smaller  amplitude  of  elastic  modes  (near  red 
dashes)  and  varies  with  elastic  parameters.  Scholte  wave  speed  increases  with  increasing  cp. 

The  left  panel  of  Fig.  3  illustrates  the  effect  of  slope  angle  on  depth-averaged  transmission  loss,  which 
serves  as  a  proxy  for  T -wave  amplitude.  The  environment  utilizes  a  Munk  deep  water  sound  speed 
profile  and  slopes  linearly  from  an  initial  ocean  depth  of  0.5  km  to  4  km  deep  over  a  range  governed  by 
the  slope  angle.  Transmission  loss  results  were  obtained  at  a  range  of  150  km.  Both  types  of  seismic 
sources  show  that  T -wave  amplitudes  increases  with  slope  angle.  In  addition,  the  transmission  loss 
results  appear  to  be  approaching  a  limiting  value  as  the  angle  increases. 

Elastic  material  parameters  impact  the  long  range  acoustic  wave  field  near  the  ocean  bottom.  For  the 
right  panel  of  Fig.  3  a  2500  m  deep  ocean  with  a  Munk  deep  water  profile  is  used.  A  10  Hz  acoustic 
source  is  placed  at  500  m  depth  and  the  signal  propagates  100  km.  A  2000  m  tall  seamount  is  located 
between  the  source  and  receiver  at  a  range  of  15  km.  When  the  acoustic  energy  interacts  with  the  elastic 
material  in  the  seamount,  both  oceanic  T -waves  and  Scholte  interface  waves  are  generated.  The  figure 
shows  wavenumber  spectra  at  zr  =  2450  m  for  increasing  values  of  cp  in  the  elastic  medium.  The  shear 
wave  speed  cs  =  0.5cp  for  all  curves.  As  cp  increases  the  amplitude  of  the  Scholte  wave  peaks  to  the 
right  of  the  ocean  wave  speed  (blue  dashes)  varies.  For  very  high  values  of  the  elastic  parameters  the 
Scholte  wave  peak  begins  to  dominate  propagating  elastic  modes.  It  is  apparent  that  the  amount  of 
energy  in  the  Scholte  wave  affects  the  amount  of  elastic  energy  that  propagates  in  elastic  body  modes 
(peaks  near  red  dashes). 
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IMPACT/APPLICATIONS 


•  Improved  modeling  capabilities  of  elastic  parabolic  equation  methods  for  underwater  acoustic 
problems  where  elastic  properties  of  the  bottom  cannot  be  ignored.  Specific  cases  include  the 
generation  and  propagation  of  oceanic  T -waves  by  seismic  sources  which  are  relevant  for 
geophysical  study  or  test  ban  treaty  monitoring. 

•  Oceanic  T -waves  and  interface  waves  are  potential  explanations  for  “deep  seafloor  arrivals”  and 
the  reception  of  acoustic  signals  in  what  may  be  otherwise  considered  a  quiet  ocean  environment 
for  monitoring. 

•  Advances  in  modeling  acoustic  propagation  in  elastic  layers  has  potential  application  in  ice 
covered  environments  where  an  elastic  layer  lies  on  top  of  the  water  column. 

RELATED  PROJECTS 

This  research  relates  to  the  separately  funded  work  of  Robert  Odom  (Applied  Physics  Laboratory, 
University  of  Washington)  regarding  the  two-way  coupled  mode  code.  It  also  relates  to  parabolic 
equation  development  by  Jon  M.  Collis  (Colorado  School  of  Mines). 
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